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Abstract 

Electromechanical systems currently offer a path to engineering quantum states 
of microwave and micromechanical modes that are of both fundamental and applied 
interest. Particularly desirable, but not yet observed, are mechanical states that 
exhibit entanglement, wherein non-classical correlations exist between distinct modes; 
squeezing, wherein the quantum uncertainty of an observable quantity is reduced 
below the standard quantum limit; and Schrodinger cats, wherein a single mode is 
cast in a quantum superposition of macroscopically distinct classical states. Also, 
while most investigations of electromechanical systems have focussed on single- or 
few-body scenarios, the many-body regime remains virtually unexplored. In such 
a regime quantum phase transitions naturally present themselves as a resource for 
quantum state generation, thereby providing a route toward entangling a large number 
of electromechanical systems in highly non-classical states. Here we show how to use 
existing superconducting circuit technology to implement a (quasi) quantum phase 
transition in an array of electromechanical systems such that entanglement, squeezing, 
and Schrodinger cats become simultaneously observable across multiple microwave 
and micromechanical oscillators. 
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Observing the quantum behaviour of macroscopic mechanical objects has received serious 
consideration in the Hterature since the turn of the century [1-8], and in 2010 O'Connell 
et al. [9] demonstrated for the first time a non-classical state of a macroscopic mechan- 
ical object by placing the dilational mode of a micromechanical resonator in a quantum 
superposition of the ground state and the first excited state. To date, however, mechanical 
entanglement, mechanical squeezing, and superpositions of " macroscopically distinct" [10] 
mechanical states (mechanical Schrodinger cats) remain unobserved. Non-local mechanical 
entanglement, first considered by Mancini et al. [11], is of fundamental interest as it would 
further elucidate the nature of the quantum-classical boundary and may provide the op- 
portunity for teleporting the quantum state of the centre of mass of a macroscopic object 
[12]. From an applied perspective, the force-sensing capability of mechanical oscillators may 
be enhanced by entanglement [13], and the expected long coherence times of mechanical 
oscillators at low temperatures may make them useful for quantum information processing 
(QIP). Squeezed mechanical states [3], where the quantum noise in one quadrature of the 
mechanical oscillator is reduced below the standard quantum limit, in addition to being of 
interest in exploring the quantuam-classical boundary can also be useful for metrology [14]. 
Finally, mechanical Schrodinger cats are relevant to tests of various models of environmen- 
tally induced decoherence [15] in mechanical resonators [16, 17] and tests of gravitational and 
non-gravitational wavefunction collapse theories (see [18] for a review or [19] for a summary). 

ECSs (see [20] for a review) are able to combine Schrodinger cats and non-local entan- 
glement in a single multipartite state such that, "The individual coherent states retain the 
desirable quasiclassical properties of coherent states but nonlocal features arise due to the 
entanglement" [21]. It is therefore highly desirable to create ECSs of mechanical resonators 
with squeezed coherent states as this would simultaneously encapsulate mechanical entangle- 
ment, squeezing, and Schrodinger cats. Bose and Agarwal proposed a scheme [22] to create 
an ECS with two nanomechanical resonators (NRs) coupled to a cooper pair box (CPB), 
and Tian and ZoUer [23] considered ECS creation in the situation of two NRs coupled to a 
trapped ion. However, both schemes are not readily scalable to larger numbers of resonators, 
do not involve squeezed states, and remain experimentally underdeveloped. 

In the photonic domain there has been much theoretical interest in ECS generation and 
use at optical frequencies [20], and the first experimental demonstration of ECS generation 
was pubhshed in 2009 by A. Ourjoumtsev et al. [25] for the case of two distant optical 
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modes. Generating ECSs in superconducting microwave circuits, however, has received httle 
attention. Superconducting circuits are an important platform for QIP and scale much more 
easily to large numbers of cavities than do optical setups. This is particularly significant for 
W-ECSs (ECSs of the form |a, 0, 0, ...0, 0) + |0, a, 0, ...0, 0) + ... + |0, 0, 0, ...0, a)) as they can 
be understood as superpositions of discrete variable W-states, for which the bipartite and 
global entanglement decay under phase damping and the global entanglement decay under 
amplitude damping have been shown to be independent of the number of qubits [26]. To 
our knowledge there have been no proposals to generate W-ECSs in microwave circuits to 
date. 

Here we present the first scheme to deterministically prepare large amplitude {\a\^ ^ 1) 
squeezed multipartite ECSs in an array of coupled superconducting microwave coplanar 
waveguides (CPWs) and to subsequently transfer them to an array of uncoupled microme- 
chanical oscillators. Mechanical mode entanglement is thereby achieved by simple state 
swapping as in the proposal of [27]. We employ the technique of entanglement generation 
via a quantum phase transition of a many-body Hamiltonian [28-30], which naturally ac- 
commodates a large number of oscillators and is particularly well-suited to the platform of 
superconducting circuits [31]. Our scheme achieves squeezing of the coherent states through 
Kerr nonlinearities in the CPWs. The resulting ECSs are able to take the following forms: 
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respectively the W-type squeezed ECS (W-sECS), W-type squeezed entangled Schrodinger 
cat state (W-sESCS), and GHZ-type squeezed ECS (GHZ-sECS), where R is the total num- 
ber of modes over which the ECS is distributed, M is the total number of CPWs, \asq) 



denotes a quadrature squeezed coherent state of amplitude the subscripts for the GHZ- 
sECS are CPW indices, and j and r are mode indices such that each mode may be a 
single fundamental mode of a CPW or the fundamental normal mode of a chain of CPWs. 
Our scheme is therefore the first to be able to distribute ECSs amongst the fundamental 
normal modes of distinct oscillator arrays [32], thereby producing non-locally superposed 
microwave/micromechanical coherent state superfluids. Also, the ECSs, in addition to being 
either purely microwave or purely micromechanical, may be distributed amongst a combina- 
tion of microwave and micromechanical oscillators. We design our scheme around existing 
technology, as we describe below, and find that parameters accessible to present-day devices 
can yield a purely mechanical version of the ECSs across ten oscillators with a multipartite 
superposition coherence of greater than sixty percent for an initial coherent state of ampli- 
tude |ap ^ 10. We also perform a numerical verification of our scheme under amplitude 
and phase damping with the master equation formalism using realistic parameters. 



SYSTEM AND MODEL 

We consider a one-dimensional chain of sites with periodic boundary conditions (Fig. 1). 
Each site consists of a superconducting microwave frequency CPW resonator, and the mi- 
crowave LC electromechanical system demonstrated by Teufel et al. in [34, 35]. A transmon 
qubit SQUID is embedded in the CPW resonator as experimentally demonstrated first by J. 
M. Fink et al. [36]. The theoretical analysis in [37] shows how this can induce a Kerr nonlin- 
earity in the fundamental microwave mode c of the CPW: He = huJcC^c — h^c^c{c^c — 1); ujc 
and respectively the frequency and Kerr constant of the mode, are dependent on the fiux 
through the SQUID loop, with cUc achieving a minimum for the fiux value that maximizes x 
and vice versa. (Higher order terms of the nonlinearity can be neglected if X ^c,o/2(c''"c), 
where c^co is the fundamental mode frequency when x = 0; see Supplementary Information). 
The electromechanical system can achieve the quantum-coherent strong-coupling interaction 
Hamiltonian Hint = —hg{o)h + afo"*"), with g as the interaction strength, which allows state 
swapping between the fundamental mechanical mode h and LC microwave mode a as recently 
demonstrated [38]. The fundamental CPW mode c is coupled to the LC mode a via the 
tunable coupler demonstrated in [39, 40], and uoc — uja when x = so that state swapping 
can occur between a and c also. The chain is formed by linking the fundamental modes Cj 
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of nearest neighbour CPWs, where j is the site index, also with the tunable coupler. All 
of the relevant dynamical time scales for our proposed system are much shorter than the 
dissipation time scales of the relevant degrees of freedom (see Fig. 1). 
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MHz - 40 MHz (~5 ns tuning time) 
0.1 MHz -50 MHz (~2 ns tuning time) 

0.1 MHz -50 MHz (~2 ns tuning time) 
MHz -1 MHz 



Cj damping rate: c 2 50 kHz 
cij damping rate: a 2 ^40 kHz 
bj damping rate: 6 2 ?^ 35 Hz 



FIG. 1. Proposed system and parameters. Grey circles denote microwave CPW (cj), microwave 
lumped LC (aj), and micromechanical (bj) oscillator modes, whose parameters are taken from refs. 
[37] and [34]. Solid black lines denote tunable inter-mode coupling. The CPWs are embedded 
with transmon SQUIDS to give a negative Kerr nonlinearity to the cj modes (ref. [37]). The 
transmon has been experimentally shown to have a Ti > 1.5 //s and a minimum of > 35 //s 
[42]. The minimum value of occurs at the flux sweet spot, but considering an asymmetric 
SQUID as in [37] allows for > 1 ms at the flux values where x is minimum and maximum 
[43]. Each CPW is tunably coupled to an LC mode of an electromechanical system of the type 
demonstrated in refs. [34, 35] and nearest neighbour CPWs by the tuneable coupler of refs. [39, 40]. 
Each electromechanical system has its own capability for tuneable coupling between its LC and 
mechanical modes [34]. (These parameters entail (c^c) < 10. A larger limit for {c^c) can be achieved 
with larger uja and maximum ujc] see Supplementary Information.) 

The model Hamiltonian for such a system, when the modes aj are all driven with equal 
and sufficient strength at frequency ujd near their red mechanical sidebands in the "resolved 
sideband" limit [? ] , can be expressed after appropriate transformations and approximations 
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where j is the site index for a total of M sites, Cm+i = Ci (periodic boundary conditions), 
5j — ujaj — ujbj — ^D,j is the detuning of the drives from the red mechanical sidebands, 
Aj — ujaj — ujcj^ 9 — 9o^/t^-) go is the bare cavity optomechanical coupling, and tid is the 
number of drive photons in the resonators of the modes aj. The Hamiltonian preserves the 
total number of quanta in the system. Because of this and the fact that the dissipation time 
scales are much longer than the relevant dynamical time scales of the unitary evolution, we 
may consider the qualitative physics of the system in an isolated and number-conserving 
picture. 

With uniform parameters {5j — 5, = A, Xj — X-) ^'j — ^")-) equation (4) is essentially 
the Hamiltonian for the one dimensional ABH with two additional degrees of freedom (cij 
and hj) linked to each lattice site. For tx! <C these additional degrees of freedom are 
effectively decoupled from the lattice, and we may make a unitary transformation to recover 
the ABH that is theoretically studied in Refs. [28, 44-52]: 



In this way our proposed system is similar to that of [53] and provides a more robust 
alternative platform to BECs in optical lattices [44, 45], coupled atom-cavity systems [30], 
and trapped ions [49] for the experimental realization of the ABH. We also note that the 
combination of quartic and quadratic terms in equation (5), given sufficient time-dependent 
control over the Xj and supplemented by coherent displacements of the fields, should 

enable one to sue this system as a universal Bosonic simulator based on a result of Braunstein 
and Lloyd [54]. In the Supplementary Information we demonstrate this for the case of 
bipartite W-ECS creation. 

In the ABH the ground state changes qualitatively as a function of the parameter r = 
—j^^z^^ where N is the total number of quanta and we assume N > 1. Any particular 
realization of the ABH can be characterized by the constants ri and T2 {t2 > ti) such that 
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for T > T2 the ground state is a superfluid, for ti > r the ground state is a Schrodinger cat- 
hke state, and the intermediate regime T2 > r > ti is transitional in nature [45]. ri ^ 0.25 
is largely independent of lattice size [45, 47], while T2 ~ [2Msiri^{7r/M)]~'^ for M > 5 and 
0.25 < r2 ^ 0.3 for 3 < M < 5 [46]. T2 = ri for M = 2 [47]. In this paper we are concerned 
with accessing the Schrodinger cat regime (r < ri), where, for decreasing values of r, the 
ground state approaches the W-state |^7v) = Tlr^j 1^)^' on-site number 

fluctuations {{cjc!jCjCj) — {c^jCj^) approach N'^^M — 1/M, and the single particle correlation 
between sites ((cjcj+i)) approaches zero [44, 45]. 

In [30] Brandao et al. suggested using an atom-cavity realization of the ABH to create a 
polaritonic (photon-atom hybrid) W-state in multiple optical cavities via adiabatic transition 
from the superfluid regime to the Schrodinger cat regime by tuning x (we note that this is not 
strictly a quantum phase transition as the thermodynamic limit for the ABH is ill-deflned 
[47]). This is possible also for the microwave quanta in our proposed system as the experi- 
mental parameter space that we specifled above allows a minimum r value of 0.0025/(A^ — 1). 
By extension, however, the linearity of quantum mechanics allows a number distribution in 
the form of a coherent state \a) = e"'^'^/^ ;^|n) (with amplitude a large enough 

such that < 0.01) to become a site-localized superposition of generalized coherent 

states: = ^ E.^ii Ur^, |0). = ^ EJi ^"'"1^/^ EZo ^'^^^H Ur^, \0)r, 

where a number-dependent phase factor e'^^'^ appears for each number state component due 
to the number-dependent frequency induced by the Kerr medium (cyclic phase shearing of 
a coherent state in a single Kerr medium was demonstrated in the recent experiment of G. 
Kirchmair et al. [56] for the case of fixed x)- For a coherent state of amplitude a localized 
on a single site, we have the exact relation (t)n^ n j ^^dt. For M > 1 the situation is com- 
plicated by the hopping term in the Hamiltonian and the fact that each number component 
sees a different value of r for the same value of x- However, our numerical investigations 
show that the different number components periodically come within a small enough phase 
difference of each other after the overall phase transition such that the coherent state at 
each site is periodically revived as a squeezed coherent state \asq) or squeezed Schrodinger 
cat state ^(e-^^/^|za,^) + e^^/^l - iasq)) to yield the W-sECS/W-sESCS/GHZ-sECS. 

A more general form of the site-localized states {"^agen) l^cf^en^^) 
= ^ X^^i l<^£^en)/ rir^j \^)r^ whcrc the lattice is divided into a total of R adjacent blocks 
with an arbitrary number of sites in each (see Fig. (2) for an example), / and r are block 
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indices, and \agen)i designates a generalized coherent state of amplitude a occupying the 
lowest energy normal mode of lattice block /: \agen)i = X^^o ^*'^''^(^!/c=o)^|0)^' 

where o\j^^^ is the creation operator for the lowest energy normal mode of block /. This cat 
state generalizes the W-ECSs studied in [57-60], which consist of coherent states confined 
to single local modes but entangled between multiple local modes, to W-ECS where any 
component of the superposition may be distributed over multiple spatially adjacent local 
modes. Such states may be thought of as non-local superpositions of superfluids. We next 
discuss the experimental protocol of how to prepare \^allT^) in our system and convert it 
into a purely mechanical state. 

PREPARATION PROTOCOL 

The experimental protocol for preparing \^all^^) in the Cj modes is illustrated in Fig. 
(2) for the case of M = 3 and i? = 2. The couplings that cross block boundaries are 
specially designated as k!" . At the beginning of the protocol, Xj = 0, and the microwave 
modes aj and Cj are all in their ground states at cryogenic temperatures and are uncoupled 
by setting k! ^ MHz. After setting k.'- = i^'Lax ^ ^0 MHz, a coherent drive is apphed 
to the lowest energy normal mode of the Cj lattice to create a coherent state of amplitude 
a, after which the drive is turned off. This results in a product coherent state of the local 
modes: I:;;;%)i|:;^)2---|:^)m- The n — 1 component of the normal mode coherent state will 
not undergo the phase transition as r is increased, so we require a sufficiently large such 
that Pi = is very small. 

The first step of the protocol adiabatically transforms the multimode coherent state of 
the Cj modes into {"^agen) by increasing x- From the discussion in the first section, we know 
that for increasing values of x/^'' the ground state with N — nis transformed at the critical 
point x//^'' 4/(n — 1) to a Schrodinger cat-like state. Considering the largest value of 
X Xmax/'^TT = 40 MHz, the bound on the smallest non-neghgible Fock state component 
l^mm) of the coherent state must be rij^in ^ ^^^^ + 1 = 5. If we impose the constraint 

Xmax 

Ylti=o K^I^)P — 0.03, we find that \a\^ > 10 is sufficient to be able to neglect number state 
components with n < 5. 

Tuning r through the transition (by changing x) to achieve 1^^^^^) involves contending 
with two dynamical processes of the system [61]. In the superfluid and intermediate regimes 

8 



1^) = -^i\(^gen)l + \agen)2 
+ |«gen)3) 

I*) = -^{\agen)l^\agen)2 
+ |agen)3) 

+|«5en)3) 

FIG. 2. Steps one, two, three, and four of the state preparation protocol ihustrated for the case of 

a total of three sites (M = 3) in two blocks {R — 2). Grey circles denote microwave CPW oscillator 

modes (cj) with negative Kerr nonlinearities of strength x- Solid black lines denote tunable intra- 

block {hi") and inter-block {i^'") couplings. The coherent state indicated in red shading can be 

site-localized (circle) or distributed over adjacent sites (oblong). Darker shading intensity reflects 

higher occupation probability. An initial coherent state is created in the lowest energy normal 

mode of the lattice, then converted into a site-localized superposition \^oigen) by adiabatically 

increasing x step one. Due to the presence of the Kerr medium, the initial coherent state \a) 

acquires a time-dependent shearing of its Wigner function and is therefore labeled as a generalized 

coherent state \agen)- In step two the blocks are decoupled by tuning n'" to zero. In step three 

X is adiabatically tuned back to zero to convert the state into the block-localized superposition 

\^al^e^^)' III step four the sites within each block are decoupled and \^all^^) is then ready to be 

transferred to the mechanical modes hj via state-swapping with modes aj. The phase shearing of 

Wgen) ccascs at the end of step three when x — ^^^d with appropriate timing the final state of 

each block corresponds to a superposition of the vacuum state |0) and either a squeezed coherent 

state \asq) or a squeezed Schrodinger cat state ^ ^\ioisq) j + e^^/^| — iasq)j). 
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(r > Ti), the hopping terms of the Hamiltonian (equation (5)) proportional to n" allow 
the quanta to redistribute into a site-localized superposition as r is tuned. Adiabaticity 
with respect to k." must therefore be maintained while x ^ i^^. From Fig. 7 of [45] 
we approximate the adiabatic condition for changing X ^ < To^^^^' which implies 
Ati ^ 0.040t2 /iS for all n > 5, where Ati is the time taken to adiabatically tune x within 
the regime r > ri. After the adiabatic tuning, x be non-adiabatically increased to 
Xmax' By increasing x while k!' is maximal, Ati is minimized and therefore the effect of the 
transmon's minimal dephasing time in the vicinity of the flux "sweet spot" can be minimized. 
On the other hand, as discussed in [44], because perfect degeneracy of the lattice sites is 
not experimentally feasible there will be non-uniformity of the lattice population induced in 
the site basis at a rate Ae/Zz, where Ae is the difference between on-site energies. This is 
relevant in the intermediate and Schrodinger cat regimes (t2 > t). If the intermediate regime 
is traversed on a time scale comparable to or longer than /i/Ae, the system will have time to 
transition to the global ground state of complete localization in the lowest energy site. We 
therefore require At^^^ <^ h/Ae^ where Atint is the time it takes to traverse the intermediate 
regime, so that the spatially uniform distribution of the quanta from the superfluid regime 
is maintained and the system ends up in a superposition of localization on each site. For 
the conservative estimate of Atint ~ Ati, this implies Ae//i < 27r x 4/t2 MHz, which can be 
satisfled for moderate lattice sizes (2 < M < 80) as reference [63] demonstrates fabrication 
capability for GHz frequency CPW resonators with variations AcUc ^ 27r x 1 MHz, while 
contributions to Ae from intersite variations in Xj can be minimized through experimental 
calibration of the control signals to keep the Xj nearly equal during step one of the state 
preparation protocol. 

Note that starting with an initial Fock state (with > 1) instead of an initial coherent 
state would yield a W-state. Our proposed scheme can therefore be an alternative to the 
one experimentally demonstrated by Wang et al. in [41], but with the added advantage of 
N being constrained by the number of sites M rather than qubit T2 times: if each site is 
tunably coupled to a qubit, M quanta can be simultaneously loaded into the array (one in 
each mode Cj via the scheme demonstrated in [62]) while = 0, after which /^'^ can be 
adiabatically tuned to Kmax to prepare the fixed number superfiuid state as the initial state 
of step one. 

After the first step, the state of the lattice is approximately {"^^agen) (assuming |ap > 10). 

10 



The third step of the protocol expands the site-locahzed \agen) of \^agen) ^^e block- 
locahzed \oLgen) of X is tuned back to (adiabatically for r > 0.25) so that the 

site-locahzed lo^^en) become block-localized superfluid |a^en) and the state of the Cj lattice 
is From the discussion of the first step we know that this must be done on a 

timescale Ata > Ati to maintain adiabaticity. Further, as discussed toward the end of the 
first section, if a is sufl&ciently large the timing may be chosen such that |a^en) ~ \(^sq) or 
Wgen) ~ ^^\'^^sq) + — ioLgq))- H thcrc is Only one block for the whole lattice 

(i? = 1), the latter case results in the GHZ-sECS of equation (3). 

In the fourth step, the intra-block couplings /^'Y27r are rapidly tuned (At4 ^ 2 ns) to 0.1 
MHz so that all the Cj modes are decoupled. This may alter the population distribution 
somewhat within each block due to imperfect diabaticity, but the superfiuid nature (inter- 
site phase coherence) of the blocks is preserved as there is no onsite interaction. 

In the final step, state transfer from Cj to aj {Atr^a ^ 0.01 /is) then aj to bj (Ats^ 0.1 /is) 
is done by pulsing and respectively. Although the mechanical modes (bj) at cryogenic 
temperatures will have significant thermal occupation, the experiment of [38] demonstrated 
that the state swap eflFectively cools the mechanical mode (bj) by transferring the thermal 
quanta to the LC mode (aj) so that the state of the mechanical mode (bj) immediately after 
the swap contains the desired state plus about one residual thermal quantum. For blocks 
with more than one oscillator, there is phase coherence between the oscillators of the block 
and the state of the block may be considered in this respect as a micromechanical superfiuid. 
We note that the conversion from microwave to micromechanical can also be site-specific 
instead of lattice- wide, resulting in an entangled microwave-micromechanical state. The 
multipartite superposition coherence after step 5b can be greater than sixty percent for 
M < 10 (see Supplementary Figure 1). 

NUMERICAL SIMULATION 

We provide a proof of principle demonstration of our proposal through numerical in- 
tegration of a master equation (see Methods) over steps one, two, and three of our state 
preparation protocol for two sites with realistic parameters as given above. We model the Ti 
and fiux dependence linearly such that Ti = 3 /iS, = 1 s when x — ^ Ti = 1.5 /iS, 

= 100 /iS when x — Xmax- A reference simulation without damping is also performed. 
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The initial state is a coherent state of amphtude \a\ — vTO in the lowest energy normal 
mode. The timing of each step respects both the minimum tuning times of each device and 
the adiabaticity constraint. Figure 3 shows the time-evolution of the fidelity {\{'^\'^w-ecs)\) 
through steps one and two, where |^^) is the damped system state and the reference state 
I^^w-Ecs) is a non-squeezed version of the W-sECS. x is tuned from zero to Xmax in 10 
ns, after which k.^^ is tuned to zero in 2 ns. The oscillations of the fidelity are indicative 
of the cyclic phase shearing of the site-localized coherent state in phase space due to the 
nonlinearity in each site of the lattice. The oscillations increase in frequency as x is increased. 

Figure 4 shows the fidelity evolution in the undamped simulation after step two if x is 
left fixed at Xmax] both the W-sECS and W-sESCS are periodically created every 27r/x = 25 
ns at the indicated separate times, and their respective Wigner functions in the insets are 
from the density matrix of the first site obtained by tracing the system density matrix over 
the second site. This indicates that by appropriately timing the tuning of x to zero in step 
three, either the W-sECS or the W-sESCS can be selected as the final state of the system: 
dt'xjnax = Jq^^ dt'x{t')^ where T* is the time at which the desired final state occurs in 
step three when x is not tuned, and T53 is the time over which x should be tuned from Xmax 
to zero in step three to achieve the desired final state. For a tuning of x linear in time we 
find that r,3 = 35 ns for the W-sECS and T.s = 10.6 ns for the W-sESCS. This is shown 
in the damped simulation results of Fig. (5) and Fig. (6), respectively. In both cases the 
fidelity oscillations decrease in frequency as x is decreased, indicating a slowing of the cyclic 
phase shearing, until the oscillations cease when x = 0. 

In Supplementary Figure 2 we show that the single particle correlator (cjc2) generally 
behaves as expected for the case of W-sECS creation, but with anomalies due to imperfect 
adiabaticity during step one, as demonstrated in Supplementary Figure 3. In Supplementary 
Figure 4, we show that less than twenty percent of the state fidelity is lost due to damping 
over steps one, two, and three of the preparation protocol. 



DETECTION 

Detection of the W-sECS/W-sESCS/GHZ-sECS may be done by swapping it into the 
aj modes and using tunably coupled superconducting qubits (as in the case of [64]) for 
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FIG. 3. Fidelity of the system state with a non-squeezed W-ECS over a simulation of steps one 
and two of the state preparation protocol for an initial coherent state of amplitude |ap = 10 in 
the lowest energy normal mode of a two site chain of CPWs. 

bipartite Wigner tomography [41] between different pairs of sites. Alternatively, the scheme 
of Tuferelli et al. [65] may be used wherein the initial state of a linear oscillator network is 
reconstructed through readout of a single qubit attached to a single oscillator in the network. 
See Supplementary Information. 



CONCLUSIONS AND OUTLOOK 

We have proposed a way to use existing superconducting circuit technology to experi- 
mentally simulate the dynamical manipulation of a many-body Hamiltonian for the purpose 
of producing a multipartite squeezed ECS of the W- or GHZ-type in the microwave domain. 
Our proposal presents a new paradigm for deterministic ECS generation that is highly rel- 
evant in the face of on-going developments in superconducting circuit-based fundamental 
quantum tests, quantum simulation, and quantum information processing. With appropri- 
ate timing of our protocol, the squeezed ECS may take the form of a non-locally superposed 
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FIG. 4. Fidelity of the system state with a non-squeezed W-ECS after step two of the preparation 
protocol when leaving x constant at Xmax- The initial state before step one is a coherent state of 
amplitude \a\'^ = 10 in the lowest energy normal mode of a two site chain of CPWs. The state of 
the system after step two periodically passes through a W-sECS (red dashed line) and W-sESCS 
(green dashed line) as revealed by the Wigner functions of the density matrix for the first site 
obtained by tracing the system density matrix over the second site at the appropriate times. 

squeezed coherent state or a non-locally superposed squeezed cat state. We have further 
shown how to combine existing electromechanical technology with the proposed supercon- 
ducting circuit to convert the multipartite photonic quantum state into a purely mechanical 
one, thereby enabling the simultaneous observation of three thus far unobserved quantum 
phenomena in mechanical systems: multipartite non-gaussian entanglement, quadrature 
squeezing below the standard quantum limit, and Schrodinger cats. This new approach to 
the preparation of entangled mechanical oscillators is one that is more readily scalable to 
many oscillators and thereby more promising for applications. Also, because our protocol 
accesses large (|ap ^ 1) squeezed coherent state mechanical Schrodinger cats, it may prove 
useful for fundamental tests in the mechanical realm. At a more general level, our proposal 
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FIG. 5. Fidelity of the system state with a non-squeezed W-ECS during step three of the prepa- 
ration protocol when x is tuned from Xmax to zero to create the W-sECS. The initial state before 
step one is a coherent state of amplitude |ap = 10 in the lowest energy normal mode of a two-site 
chain of CPWs. The timing is chosen such that the stationary state of the system after the tuning 
period is a W-sECS, as indicated by the Wigner function of the first site obtained by tracing the 
system density matrix over the second site at the end of the tuning period. The fidelity oscillations 
decrease in frequency as x is decreased and cease when x — 0- 

highlights the use of quantum phase transitions as a resource for fundamental and applied 
studies in engineered quantum systems, especially in the case of multipartite entanglement. 



METHODS 

The numerical simulations are performed by using the fourth order Runge-Kutta inte- 
gration method on the master equation: 

i i ^ 
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FIG. 6. Fidelity of the system state with a non-squeezed W-ECS during step three of the prepara- 
tion protocol when x is tuned from Xmax to zero to create the W-sESCS. The initial state before 
step one is a coherent state of amplitude |ap = 10 in the lowest energy normal mode of a two-site 
chain of CPWs. The timing is chosen such that the stationary state of the system after the tuning 
period is a W-sESCS, as indicated by the Wigner function of the first site obtained by tracing the 
system density matrix over the second site at the end of the tuning period. The fidelity oscillations 
decrease in frequency as x is decreased and cease when x — 0- 



where p is the density matrix for the CPW chain, H is the Hamiltonian of the ABH model 
from equation (5), T>[cj]p = Cjpc^- — c^-Cjp/2 — pc^-Cj/2 is the amplitude damping operator, 
Q[cj\p = c^jCjpc^jCj — {c^-CjYp/2 — p(cjcj)^/2 is the phase damping operator, 1/Ti is the 
amplitude damping rate, and l/T^^ is the phase damping rate. The Hilbert space at each 
site is truncated at rimax = 20, and a timestep of size 10~^^ s is used. 
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I. CONSTRAINT ON COHERENT STATE AMPLITUDE TO NEGLECT 
HIGHER ORDER TERMS OF CPW NONLINEARITY 



The full expression for the nonlinearity (Hnl) introduced into the CPWs by the transmon 
qubits is given in equation 27 of [1]. Looking only at the fundamental mode, and using the 
relations $o = h/2e, Ej = {^f/Lj, Lj = L'/rn, and L' = 2Ey {001^6^) we find 

*- = E^(^)"'»'^.o[<='+^-. (1) 

where ^/ (0 < 77/ < 1) is the inductive participation ratio dependent upon the flux through 
the SQUID loop and a;c,o is the fundamental mode frequency in the absence of the transmon. 
After the rotating wave approximation, to be able to neglect nonlinear terms higher than 
{c^cy we require 



or 




where we use the fact that greater than 95% of a coherent state with amplitude a occupies 
number states with n < 1.5|ap for lap > 10. Using x = 40 MHz, cj^o — 7.5 GHz, and 
lap = 10, we find that the constraint is satisfied. We note that larger values of x and/or 
jap can be accommodated by increasing the bare CPW fundamental mode frequency cj^o 
as long as the lumped LC oscillator frequency cUa matches lUc^q. Larger values of x in turn 
allow for larger values of /^^^x therefore smaller Ati and Ats, thereby increasing the 
quantum coherence of the state when transferred to the mechanical modes. 



II. BOSONIC SIMULATOR DEMONSTRATION FOR BIPARTITE W-ECS CRE- 
ATION 

Consider for the case of two sites the Hamiltonian given in equation (5) of the main 
article, 

f = -^clcMci - 1) - ^4C2(4C2 - 1) - K"{t){clc2 + Ci4), (4) 

with the explicit assumption that each of the quartic non-linear terms and the quadratic 
coupling term may be subject to independent control. 
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We begin by assuming that the Kerr nonhnearity is set to zero. Then both normal 
modes, Ai = (ai + = (ai — a^j^^ of the coupled cavity system can each be 

coherently driven to a steady-state coherent state. The state of the local modes is likewise 
a product coherent state, and by a suitable choice of the normal mode driving fields it 
becomes \oi)ax \ — ^a^- the first step we set ^' to zero, to decouple the cavities, and 
the Kerr nonlinear term for cavity a\ is to a non-zero value for a fixed time T\ such that 
= TT. The Kerr nonhnearity for mode remains set at zero. This transforms the local 
modes as 

|^^)ai| - (^)a2 ^ [e~'^'%a)a^ + e'^l^\ - ia)a,) I - a)a, (5) 

In the next step we turn on the coupling between the cavities for a fixed amount of time T2 
so that /^''T2 = 7r/4, the resulting state is 

e-^^/4|zV2a)„jO),, + e'^/^\0)a,\ - V2a)a, (6) 

While this is not precisely a W-ECS it is a non-classical entangled cat state. Using suitable 
phase rotations and displacements however it can easily be transformed into a W-ECS state. 

III. REMAINING COHERENCE FRACTION AFTER STATE PREPARATION 

We can estimate the quantum coherence of the ECS after it is prepared and transferred to 
the mechanical modes. Using the formula 7^ = ^ + (see [2]), an average Ti ^ 2 /is, and 
an average dephasing time of T^^ 0.1 s over the relevant fiux range, we estimate the average 
decoherence over steps one, two, three, and four for a coherent superposition |n)|0) + |0)|n) 
across two of the Cj modes to be D^^^M^) = e-(^*i+^*2+^*3+At4)(n/(2xio-6)+nV(iooxio-6)^ 

The decoherence D^aip) — 6"^*^"/^°^^ over Ats^ entails an average of the decay rates of the 
CPW and the LC oscillator (we assume dephasing to be negligible for both modes when 
^ = 0): Ton = 2/n(7c + 7a) 5 = 3.42/n jis. The decoherence ^^(n) = ^-^t^h/Ton ^^^^ 
At55 entails an average of the decay rates of the LC oscillator and the mechanical mode (we 
assume dephasing to be negligible for these modes): T^^ = "^/^{penvlh + la) s = 7.67/n 
where Uenv — 50 is due to the mechanical mode thermal environment at 25 mK [3]. In 
Fig. (1) we plot the total remaining quantum coherence of a bipartite W-ECS formed in 
a two site lattice with coherent state amplitude a after all of the preparation stages under 
the approximation that the coherent state is already in ECS form at the beginning of step 
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one: Zl^^o K^I^)P^i,2,3,4(^)^5a(^)^56(^) vs. \a\. We use Ati = .02 /xs, At2 = .002 /xs, 
Ata = .02 /xs, At4 = .002 /is, Ats^ = .01 /xs, and Ats^ = .16 /xs. The values of Ati and 
At2 accommodate M < 10. The plot indicates that our proposed system and protocol can 
achieve at each site a purely mechanical cat state with greater than sixty percent quantum 
coherence and an average centre of mass spatial separation of a few times the mechanical 
zero point motion. 




2 4 6 8 10 



\a\ 



FIG. 1: Approximate remaining fraction of quantum coherence as a function of normal mode 
coherent state amplitude \a\ of initial M-partite W-ECS (M ^10) after transfer to mechanical 
modes. The calculation is done under the approximation that the coherent state is already in ECS 
form at the beginning of step one. 



IV. BEHAVIOUR OF (cjcs) FOR THE CASE OF W-SECS CREATION 

The single particle correlator (c|c2) between sites one and two is plotted in Fig. 2 over 
all three steps for the case of W-sECS creation plus an additional 5 ns free evolution period 
after step three. The steady overall decrease of (c|c2) indicates increasing site-localization 
of the initial coherent state (the perfect ECS has complete site-localization and a single 
particle correlator value of zero). The intermittent oscillations appear to have a frequency 
proportional to x, which indicates that they are a signature of dynamics between the different 
eigenmodes of the lattice. This arises because the nonlinear term of the Hamiltonian induces 
scattering between the normal modes (as easily seen by Fourier transforming the local mode 
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operators Cj), and because the tuning of x step one is not perfectly adiabatic. This is 
verified below. The W-type distribution of the coherent states is not affected by this because 
all of the low energy states in the regime r < 0.25 are of the W-state form with different 
local phase variations only [4]. 





(4C2) 




t^" 2 




X 2 




45 
40 
35 
30 
25 
20 
15 
10 
5 



0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 



FIG. 2: (c|c2) over steps one, two, three, and a 5 ns period after step 3. The oscillation pulses are 
discussed below. The ideal ECS has a value of zero, indicating complete site-localization. 

To show that the amplitude of the intermittent oscillations decreases as Ati is increased, 
we perform three undamped numerical simulations of steps one and two of the state prepa- 
ration protocol followed by a period of unitary evolution with x = Xmax- Each simulation 
uses different values of Ati in order to test the dependence of the oscillations of {c\c2) on 
the adiabaticity of step one (where x is tuned from zero to Xmax) and the results are shown 
in Fig. 3. The oscillation amplitudes clearly decay with increasing Ati, indicating that the 
oscillations are due to the imperfect adiabaticity of step one. 
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FIG. 3: Undamped simulation of (c|c2) over steps one and two, followed by a period of unitary 
evolution with x — Xmax- si, Ati = 0.02 fis. b, Ati = 0.04 fis. c, Ati = 0.08 /xs. With increasing 
Ati the adiabaticity of step one is improved and the oscillations are consequently increasingly 
damped. The steep drop after Ati is due to step two, where the sites are decoupled by tuning k,^^ 
to zero. Other than Ati, the simulation parameters are the same as for the undamped simulation 
of Fig. 5 in the main text. 

6 



V. EFFECT OF AMPLITUDE AND PHASE DAMPING ON THE FIDELITY OF 
THE CREATED STATE 

To investigate the effect of the amphtude and phase damping on the state preparation, 
we consider the fidehty F{p^p') = [tr-^y/pp'y/p]^ of the damped system state (p) with the 
undamped system state (p^) through aU steps of the preparation protocol for the W-sECS. 
As F(p, p') is computationaUy very expensive to compute, we calculate instead the upper 
bound, given by F(p, p^) < trpp' + ^(1 — trp'^){l — trp''^)^ and the lower bound, given by 
F{p,p') > trpp' + V2^{trpp'y - trp^p^^ [5]. Both are shown in Fig. 4, which clearly 
indicates that less than twenty percent of the fidelity is lost due to damping over steps one, 
two, and three of the preparation protocol. 
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FIG. 4: Upper bound (dashed purple) and lower bound (solid black) of the fidelity of the system 
state p with the undamped reference system state p' over the preparation of the W-sECS. 
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VI. STATE DETECTION VIA A SINGLE QUBIT 



After creation of the mechanical/microwave W-sECS/W-sESCS/GHZ-sECS in our sys- 
tem, all couplings k! ^ k!' ^ and g could be turned on while leaving x = to make the 
network linear. In this case, the scheme of Tuferelli et al. [6] would allow for a single qubit 
tunably coupled to any microwave mode in our system to reconstruct the initial mechani- 
cal/microwave W-sECS/W-sESCS/GHZ-sECS of the system. One caveat of this approach 
is that it requires Markovian master equation modelling of the dissipative environment for 
all of the oscillators of the network. Because the correct such model for the mechanical 
dissipation is unknown, the bipartite Wigner tomography method may need to be used to 
first verify which dissipative model to use before applying the method of Tuferelli et al. to 
read out the entire network. 
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